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Abstract 

Physical systems, characterized by an ensemble of interacting elementary constituents, can be represented 
and studied by different algebras of observables or operators. For example, a fully polarized electronic sys- 
tem can be investigated by means of the algebra generated by the usual fermionic creation and annihilation 
operators, or by using the algebra of Pauli (spin- 1/2) operators. The correspondence between the two al- 
gebras is given by the Jordan- Wigner isomorphism. As we previously noted similar one-to-one mappings 
enable one to represent any physical system in a quantum computer. In this paper we evolve and exploit 
this fundamental concept in quantum information processing to simulate generic physical phenomena by 
quantum networks. We give quantum circuits useful for the efficient evaluation of the physical properties 
(e.g, spectrum of observables or relevant correlation functions) of an arbitrary system with Hamiltonian H. 

PACS numbers: Pacs Numbers: 3.67. Lx, 5.30.-d, 
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I. INTRODUCTION 



A fundamental concept in quantum information processing is the connection of a quantum 
computational model to a physical system by transformations of closed operator algebras. The 
concept is a necessary one because in quantum mechanics each physical system is naturally as- 
sociated with a language of operators (for example, quantum spin- 1/2 operators) and thus to an 
algebra realizing this language (e.g., the Pauli spin algebra generated by a family of commuting 
quantum spin- 1/2 operators). Any quantum system defined by an algebra of operators generated 
by a set of "basic" operators can be considered as a possible model of quantum computation flip. 
The remarkable fact is that an arbitrary physical system is simulatable by another physical system 
(or quantum computer) whenever isomorphic mappings (embeddings) between the two operator 
algebras exists. In each such case, an important problem is to determine whether the simulation 
is efficient (polynomial resource overhead) in terms of the "basic" generators. For example, a nu- 
clear spin (NMR) quantum computer is modeled as a collection of quantum spin- 1/2 objects and 
described by the Pauli algebra. It can simulate a system of 4 He atoms (with space discretized by a 
lattice) represented by the hard-core bosonic algebra, and vice versa. In this case, the simulation is 
efficient. Figure [1] summarizes this fundamental concept by giving a variety of proposed physical 
models for quantum computers and associated usable operator algebras. If one of these systems 
suffices as the universal model of quantum computing, the mappings between the operator alge- 
bras establish the equivalence of the other physical models to it. This is one's intuitive expectation, 
and has a well-established mathematical basis [|3p. 

The mappings between algebras, between an algebra and a physical system, and between phys- 
ical systems are necessary in order to be able to simulate physical systems using a quantum com- 
puter fabricated on the basis of another system. However, this does not imply that the simulation is 
efficiently implementable. As we have previously discussed [[I]], efficient quantum computation in- 
volves more than having the ability to represent 2^ different items of classical information so that 
the algebra of N quantum bits (qubits) can be isomorphically represented and quantum parallelism 
can be exploited. It is also insufficient for the mapping between operator algebras to be easily and 
perhaps efficiently formalized symbolically. For example, the physical system consisting of one 
boson in 2^ modes is described using the language of "transition" operators that move the boson 
from one mode to the other. Formally, the Pauli matrices on N qubits can be easily represented 
using the transition operators, but the one-boson system is no more powerful than classical wave 
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mechanics. This means that unless quantum computers are not as powerful as is believed, there is 
no efficient simulation of qubits by the one-boson system. 

To be useful as a physics simulation device, a quantum computer must answer questions about 
physical properties associated with real physical systems. These questions are often concerned 
with the expectation values of specific measurements of a quantum state evolved from a specific 
initial state. Consequently, the initialization, evolution, and measurement processes must all be 
implementable with polynomial scaling [jlj]. Often it is difficult to do. Further, some classes of 
measurements, such as thermodynamic ones, still lack well-defined workable algorithms [Q]. 

On a classical computer, many quantum systems are simulated by the Monte Carlo method [§]. 
For fermions, the operation counts of these Monte Carlo algorithms scale polynomially with the 
complexity of the system as measured by the number of degrees of freedom, but the statistical 
error scales exponentially (in time and in number of degrees of freedom), making the simulation 
ineffective for large systems. A quantum computer allows for the efficient simulation of some 
systems that are impractical on a classical computer. In our recent paper [[]]] we discussed how 
to simulate a system of spinless fermions by the standard model of a quantum computer, that is, 
the model expressed in the language and algebra of quantum spin- 1/2 objects (Pauli algebra). We 
also discussed how to make certain physically interesting measurements. We demonstrated that 
the mapping between algebras is a step of polynomial complexity and gave procedures for initial 
state preparation, evolution, and certain measurements that scaled polynomially with complexity. 
The main focus of the paper however was demonstrating that a particular problem for simulating 
fermions on a classical computer, called the dynamical sign problem, does not exist on a quantum 
computer. We are aware of at least one case where the sign problem can be mapped onto an NP- 
complete problem [§]. This is the 3-SAT problem [|7|]. Therefore, one cannot yet claim that a 
quantum computer can solve "all" sign problems, otherwise one would claim that one is solving 
all NP-complete problems and this has not been rigorously established. 

In this paper we continue to explore additional issues associated with efficient and effective 
simulations of physical systems on a quantum computer, issues which are independent of the par- 
ticular experimental realization of the quantum computer. We seek to construct quantum network 
models of such computations. Such networks are sets of elementary quantum gates to which we 
map our physical system. For simplicity, we discuss these issues relative to simulating a system 
of spin- 1/2 fermions by the standard model of quantum computing. Our discussion has obvi- 
ous applications to the simulation of a system of bosons (or any other particle statistics or, in 
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mathematical terms, any other operator algebra). Specifically we address issues discovered in our 
attempt to implement a (classical) simulator of a network-based quantum computer and to conduct 
a quantum computation on a physical system (NMR) with a small number of qubits. On a classical 
computer the number of qubits simulatable is limited by the exponential growth of the memory 
requirements. Physically, we can only process information experimentally with systems of a few 
qubits. Having the simulator permits a comparison between theory and experiments likely to be re- 
alizable in the near future. Overall, the main problems we address are how to reduce the number of 
qubits and quantum logic gates needed for the simulation of a particular physical phenomena, and 
how to increase the amount of physical information measurable by designing efficient quantum 
algorithms. 

We organized the paper in the following manner: In Section II we summarize the quantum 
network representation of the standard model of quantum computation, discussing both one- and 
multi-qubit circuits. Then we summarize the connection between the spin and fermion represen- 
tations. In Section we first discuss the initialization, evolution, and measurement processes. 
In each case we define procedures simpler than the ones presented in our previous paper, greatly 
improving the efficiency with which they can be done. Greatly expanded are the types of measure- 
ments now possible. For example, besides certain correlation functions, the spectrum of operators, 
including the energy operator, can now be obtained. Our application of this technology to a system 
of fermions on a lattice and the construction of a simulator is discussed in Section [IV]. The Hub- 
bard model is used as an example. We conclude with a summary and a discussion of areas needing 
additional work. The appendices contain technical points about the preparation of coherent and 
correlated states and the use of the discrete classical Fourier transformation. 

II. QUANTUM NETWORK REPRESENTATION OF PHYSICAL PHENOMENA 

It is the formal connection between models of computation and physical systems described 
in the Introduction that allows one to simulate quantum phenomena with a quantum computer. 
Simulation is realized through a quantum algorithm that consists of unitary operations and mea- 
surements. One of the objectives is to accomplish simulation efficiently, i.e, with polynomial 
complexity. The hope is that quantum simulation is "more" efficient (less resources) than classical 
simulation and there are examples that support such hope [0. In the following subsections we 
summarize the main concepts in the representation of physical phenomena by quantum networks. 
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A. Standard Model 



In the standard model of quantum computation, the quantum bit, or qubit, is the fundamental 
unit. A qubit's state |a) is a linear combination of the states |0) and |1) (e.g, a spin 1/2 with 

|0) = |T), |l> = li»: 

|a)=a|0) + 6|l>, (1) 

where the complex numbers a and b are normalized to unity: |a| 2 + |5| 2 = l. 

Assigned to each qubit are the identity matrix 11 and the Pauli matrices a x , a y and a z : 

1= (! !)'"" = C o)"" = G -i)' <2> 

or equivalently 11, a± = ^(a x ±ia y ), and a z . In this particular representation, the states |0) and |1) 
are the vectors: 

|0)= (]j and|l)= (°) , (3) 

and the Bloch-Sphere (Fig. [|) provides a convenient three-dimensional real space representation 
of the single qubit state |a), which can be parametrized as |a) = cos §|0) + e lLp sin 

For a system of n qubits, the mathematical representation of the standard model is defined by a 
closed *-algebra (Pauli algebra) generated by the operators aj^ (/i = x, y, or z) that act on the j th 
qubit: 

n factors 

A, 



o£ = 1 ® 1 ® • • • ® cr^ ® • • • <g) 1 , 

j th factor 

where ® represents a Kronecker product. From these definitions, the resulting commutation rela- 
tions are 

= 25, v (4) 
= ^jke^xai , (5) 

where [A, _B]± = AB±BA, and e Ml/ A is the totally anti- symmetric Levi-Civita symbol. The time 
evolution of an n qubit system is described by the unitary operator U (t) = e~ tHt , where H repre- 
sents the time-independent Hamiltonian of the system. In turn, U(t) is easily expressible in terms 
of the Pauli matrices cr^ since they and their products form an operator basis of the algebra. 
The most general unitary operator U on a single qubit can be written as 

U = e ia R z (P)Ry(l)Rz(S), (6) 



where a, (5, 7, and 5 are real numbers, and R^) = e _t 2 <jM are rotations in spin space by an angle 
•d about the \i axis. Although this decomposition is not unique, it is important because any one 
qubit evolution is seen to be a combination of simple rotations (up to a phase) about the fi = x,y 
or z axis. 

In multi-qubit operations, any unitary operation U can be decomposed (up to a phase) as U = 
~[ Ui, where Ui are either single qubit rotations R^fl) in the n-qubit space or two qubit interactions 

R z i,z k {u) = e tu)<Jza * in the same space (u is a real number) [fjj, §]. These one qubit rotations and 
two qubit interactions constitute the elementary gates of the quantum computer in the network 
model. 



B. Quantum Network 

We now describe some common one and two qubit gates, some quantum circuits, and one 
pictorial way to represent them. The motivation for this elementary subsection is to prepare the 



grounds for the quantum network simulation of a physical system developed in Section [III] which 
is more technically involved. 

The goal is to represent any unitary operation (evolution) as a product of one and two qubit 
operations. Although here we use the algebra of the Pauli matrices (standard model), for a different 
model of computation we should change the set of elementary gates, but the general methodology 
remains the same. For instance, if the evolution U(t) = e~ tHt is due to the Hamiltonian 

H = H x + H y = a a\a 2 z ■ ■ ■ a{~ l ai + f) a\a 2 z ■ ■ ■ a{- l a j y , (7) 

where a and j3 are real numbers, we write U(t) as e - lB xt e -iH v t b ecause [H x ,H y ]_ = 0. To 
decompose this into one and two qubit operations, we take the following steps: We first note that 
the unitary operator 

t/i = e^ = -J=[i + ^J] (8) 
takes o\ — > o\, i.e., \j\o\XJ\ = a\, so U\e %ol<Jz U\ = e lo " Tx . Next we note that the operator 

takes a\ —> o\a 2 z , so £7 2 t e ia<T ^£y 2 = e mf7 ^. Then we note that 

U 3 = 



takes o-y(j 2 z — > —ala 2 a z . By successively similar steps we easily build the required string of 
operators: o\o 2 z ■ ■ ■ o{~ x al. and also e iaa i cr l- a i 1<J ^ ( U p to a global phase): 

Ul • • • Ulule^UxUi ■■■U k = e ™*° a .~°i~ 1 <d (9) 

where the integer k scales polynomially with j (in this particular case the scaling is linear). In a 
similar way, we decompose the evolution e~ lHyt . Multiplying both decompositions, we have the 
total decomposition of the evolution operator U (t). See QTO , [TT| ] for complete treatments of these 
techniques. 



1. Single Qubit Circuits 

In Fig. |3]a we show examples of several elementary one qubit gates. (Notice that e ie °~^ = 
R^—26).) Each gate applies one or more unitary operations R^fi) to the qubit (the gates 
apply a R^ir) rotation up to a phase: = ie~ l ^ CTM ). Also, in Fig. we show the Hadamard gate 
H. The action of this gate on the state of one qubit |a) is: 



H 



|0)-|+) = (|0) + |l))/v/2 
|l)-h) = (|0)-|l))/v / 2 

In this way, the Hadamard gate admits the matrix representation: 



h ^G-i'' <10) 

In terms of the Pauli matrices 

H = -j=[a x + <r g ] = ie-^ ax e~^ ay . (11) 

In Fig. |]a we show the decomposition of the H gate into single qubit rotations, and its application 
to the Bloch-Sphere representation of the state |+) is shown in Fig. |]b. The convention for quan- 
tum circuits is each horizontal line represents the time evolution of a single qubit and the time axis 
of the evolution increases from left to right. 



2. Multiple Qubit Circuits 

We now give examples of multi-qubit operations. Again the goal is to represent them as a 
combination (up to a phase) of single qubit rotations -R M ($) and two qubit interactions R z j jZ k (uj) = 

7 



e' 



Jala 



(the gate for the R z j z k (u) is shown in Fig. |3]b). To illustrate this, we consider the circuit 



shown in Fig. |[ This is a two qubit controlled-NOT (C-NOT) gate which acts as follows: 



C-NOT < 



|oo>- 


-|oo) 


|01>- 


- |oi> 


|10>- 


-111) 


Jll>- 


+ |10) 



Here, the first qubit is the control qubit (the controlled operation on its state |1) is represented by 
a solid circle in Fig. |5|). We see that if the state of the first qubit is |0) nothing happens, but if the 
first qubit is in |1), then the state of the second qubit is flipped. Because a\ is the unitary operator 
that flips the second qubit (see Fig. |J), the decomposition of the C-NOT operation into one and 
two qubits interaction is 

C-NOT: e i %e~ i % a *e~ i * a *e i % a *' J * = e i ?e~ i %°'e~ i %°*e i %tfe~ i ?'^' 7 *e~ i % a % . (12) 

From Eq. [T2| we can see that a single controlled operation becomes a greater number (in this 
case 4) of one and two qubits operations. In Fig. ^] we also show the circuit representing this 
decomposition, while in Fig. |^ we show the C-NOT gate applied to the state 1 10) in the Bloch- 
Sphere representation. Because of the control qubit being in the state |1), the second qubit is 
flipped. 

A generalization of the C-NOT gate is the controlled-t/ (C-U) gate, where U is a unitary oper- 
ator acting on a multi-qubit state |W S ): 



C-U 



|0) a ® |* a > -> |0) a <8> \V S ) 

|l)a® |*,) -> |l) a ® \U\Vs) 



Mathematically, for U(t) = e lQt (Q is Hermitian), the operational representation of the C-U gate 
is: U(t/2)U(t/2)- a * (U(t)-°** = e^®^ 1 ), where a is the control qubit (Fig. 0a). Similarly, one 
can use |0) a as the control state to define the C-U' gate illustrated in Fig. [7]b. In order to describe 
the C-U and C-U' gates as a combination of single qubit rotations and two qubits interactions, 
we have to decompose the operators U(t/2) and U(t/2) a ** into such operations. C-U can then 
be expressed as a sequence of conditional one and two qubit rotations. The latter can be further 
decomposed into one and two qubit rotations using the techniques of [Bp. 
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C. Spin-Fermion connection 



To simulate fermionic systems with a quantum computer that uses the Pauli algebra, we first 
map the fermionic system into the standard model [[I], [T2p . The commutation relations for (spinless) 
fermionic operators aj and at (the destruction and creation operators for mode j) are 

[a],a k } + = S jk , (13) 
[a},4] + = 0. (14) 

We map this set of operators to one expressed in terms of the a^'s in the following way: 

a i -> (n -°i J °- = ■ ■ ■ ^rv: , 
°5 -> (n < = i-^ 1 ^ 2 • ■ ■ < x < ■ 

Obviously, for the fermionic commutation relations to remain satisfied, the operators cr-? must 
satisfy the commutation relations of the Pauli matrices, so a representation for the operators <ji are 
the Pauli matrices. 

The mapping just described (indeed it induces an isomorphism of *-algebras) is the famous 
Jordan-Wigner transformation [O]. Using this transformation, we can describe any fermionic uni- 



tary evolution in terms of spin operators and therefore simulate fermionic systems by a quantum 
computer. Although the mapping as given is for spinless fermions and for one-dimensional sys- 
tems, it extends to higher spatial dimensions and to spin- 1/2 fermions by re-mapping each "mode" 
label into a new label corresponding to "modes" in a one-dimensional chain. In other words, if we 
want to simulate spin- 1/2 fermions in a finite N x x N y two-dimensional lattice, we map the label 
of the two-dimensional lattice to an integer number S, running from 1 to 2(N X x N y ). S identifies 
a mode in the new chain: 

- a 5 - -*i ) - s = (-i)*- 1 ■ ■ ■ > 

<W - 4 - (it 4 = ° X A ■ ■ • *fM , (is) 

where the d(j,k)\a an d a \jk)-a are *- ne fermionic spin- 1/2 operators in the two-dimensional lattice 
for the mode (j, k) and for z-component of the spin a (a — ±|), and a s and ajy are the spinless 



fermionic operators in the new chain. In our case, the modes are the sites and the label (j, k) 
identifies the X-Y position of this site (j, k e [1, N Xjy \). The label (j, k); a maps into the label S 
(Fig. |J) via 

S = j + (k - 1)N X + (i - a)N x N y , (16) 

This is not the only possible mapping to a two-dimensional lattice using Pauli matrices [Q, |14|,|15[], 
but it is very convenient for our simulation purposes. 

III. QUANTUM NETWORK SIMULATION OF A PHYSICAL SYSTEM 

Like the simulation of a physical system on a classical computer, the simulation of a physical 
system on a quantum computer has three basic steps: the preparation of an initial state, the evolu- 
tion of the initial state, and the measurement of the physical properties of the evolved state. We will 
consider each process in turn, but first we note that on a quantum computer there is another impor- 
tant consideration, namely, the relationship of the operator algebra natural to the physical system 
to the algebra of the quantum network. Fortunately, the mappings (i.e., isomorphisms) between 
arbitrary representations of Lie algebras are now known []3|]. Section pTC] is just one example. To 
emphasize this point, the context of our discussion of the three steps will be the simulation of a 
system of spinless fermions by the standard model, which is representable physically as a system 
of quantum spin 1/2 objects. 

A. Preparation of the Initial State 

The preparation of the initial state is important because the properties we want to measure 
(correlation functions, energy spectra, etc.) depend on it. As previously discussed [[[]], there is a 
way to prepare a fermionic initial state of a system with iV e spinless fermions and n single particle 
modes j, created by the operators a] ( creation of a fermion in the mode j). In the most general 
case, the initial state is a linear combination of Slater determinants 

Ne 

i $ «)=ri 6 ]i vac )' < iv) 

described by the fermionic operators bj and fet which are related to the operators aj and a] via a 
canonical (unitary) transformation. Here |vac) is the vacuum state (zero particle state). To prepare 
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\& a ) one can look for unitary transformations U m such that 



Ne 



\$ a ) = ^ H U m |vac) , (18) 

m=l 

where 7 is a phase factor. To perform these operations in the standard model we must express the 
U m in terms of Pauli matrices using the Jordan- Wigner transformation. (We can do the mapping 
between the Pauli operators and the a, operators or between the Pauli operators and the bj opera- 
tors. In the following we will assume the first mapping since this will simplify the evolution step.) 
One can choose U m = e~ lHmt such that H m is linear in the b rn and operators []Tp. We have to 
decompose the U m into single qubit rotations and two qubit interactions -R M ($) and R z j^ z k(uj). To 
do this, we first decompose the U m into a products of operators linear in the b m or 6^; however, 
this decomposition does not conserve the number of particles. The situation appears complex. 



Simplification occurs, however, by recalling the Thouless's theorem QI6Q which says that if 

N e 

|0) = J]4l VaC ) < 19 ) 



and M is a n x n Hermitian matrix, then 

Acft Ma 



N e 

n&J|vac>, (20) 



where a) = (a\, ■ ■ ■ , a^) and 



P = e lM a f . (21) 



From Eq. ^T]the operator e lM (formally acting on the vector of aj's) realizes the canonical trans- 
formation between cij and bj. 

Thouless's theorem generalizes to quantum spin systems via the Jordan- Wigner transformation. 
This theorem allows the preparation of an initial state by simply applying the unitary operator 
e ia Ma t0 a "b oot U p" state polarized with each qubit being in the state |0) or |1). Indeed, for an 
arbitrary Lie operator algebra the general states prepared in this fashion are known as Perelomov- 
Gilmore coherent states JT7| ] . 

The advantage of this theorem for preparing the initial state instead of the method previously 
described [jl]] is that the decomposition of the unitary operator e mtA/a can be done in steps, each 
using combinations of operators aja\ and, therefore, conserving the number of particles. Once 
the decomposition is done, we then write each operator in terms of the Pauli operators to build a 
quantum circuit in the standard model. (See Appendix A for a simple example.) 
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A single Stater determinant is a state of independent particles. That is, from the particle per- 
spective, it is unentangled. Generically, solutions to interacting many-body problems are entangled 
(correlated) states, that is, a linear combination of many Slater determinants not expressible as a 
single Slater determinant. In particular, this is the case if the interactions are strong at short ranges. 
In quantum many-body physics, considerable experience and interest exists in developing simple 
approaches for generating several specific classes of correlated wave functions \$Wj\. In Appendix 
A we illustrate procedures and recipes to prepare one such class of correlated (entangled) states, 



the so-called Jastrow states [16]. 



B. Evolution of Initial State 

The evolution of a quantum state is the second step in the realization of a quantum circuit. The 

goal is to decompose this evolution into the "elementary gates" R^fi) and R z j^ z k(u). To do this 

for a time-independent Hamiltonian, we can write the evolution operator as U (t) = e~ lHt , where 

H = ^ Hi is a sum of individual Hamiltonians Hi. If the commutation relations [Hi, Hi]_ = 
i 

hold for all / and I', then 

U(t)=Hu l (t) = l[e- iH > t . (22) 

i i 

In this way, we can then decompose each Ui(t) in terms of one and two qubits interactions, using 
the method described in Section PB| . 

In general, the Hamiltonians Hi for different / do not commute and the relation Eq. cannot 



be used. Although we can in principle exactly decompose the operator U (t) into one and two 
qubit interactions ffl, ^p, such a decomposition is usually very difficult. To avoid this problem, 
we decompose the evolution U(t) = JJ e~ lHAt using the the first-order Trotter approximation 
(t = MAt): 

U(At) = e~ iHAt = e ~ l ^ HlAt = Y[ e~ lH ' At + 0((At) 2 ) . (23) 

i 

Then, for At — > 0, we can approximate the short-time evolution by: U(At) ~ Y\e~ lHlAt . In 

i 

general, each factor is easily written as one and two qubit operations (Section PB[ ). 

The disadvantage of this method is that approximating the operator U(t) with high accuracy 
might require At to be very small so the number of steps e~ lHlAt and hence the number of quantum 
gates required becomes very large. To mitigate this problem, we can use a higher-order Trotter 
decomposition. For example, if H = K + V, we then use the second-order Trotter approximation 
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to decompose the evolution as U(t) = Yl e lHAt with (second-order decomposition) 

3 

e -iHM = e - iK % e -iVM e -iKM + 0((At) 3 } ? (24) 

= e -* v ^e- iKAt e- tV ¥ + 0((At) 3 ) . (25) 



Other higher-order decompositions are available Ql 8Q. 



C. Measurement of Physical Quantities 

1. One-Ancilla Qubit Measurement Processes 

The last step is the measurement of the physical properties of the system that we want to study. 
Often we are interested in measurements of the form (WV), where U and V are unitary operators 
[|T]] . We refer to Ref . [|T]] for a description of the type of correlation functions that are related to 
these measurements. See also [ JT9| ] for an application and variation of these techniques. Here, we 
simply give a brief description of how to perform such measurements. 

First, we prepare the system in the initial state |^ ) an d adjoin to it one ancilla (auxiliary) 
qubit a, in the state |+) = (|0) + |l))/\/2. This is done by applying the unitary Hadamard 
gate to the state |0) (Fig. |]). Next, we make two controlled unitary evolutions using the C-U 
and C-U' gates. The first operation V evolves the system by V if the ancilla is in the state |1): 
V = 1 0) (0 1 ® i + 1 1 ) ( 1 1 ® "K. The second one U evolves the system by U if the ancilla state is 1 0) : 
U = |0}(0| ® U + ® 1. (V and U commute.) Once these evolutions are done, the expectation 
value of 2o\ = <j\ + ia^ gives the desired result. This quantum circuit is shown in Fig. |9|. Note 
that the probabilistic nature of quantum measurements implies that the desired expectation value 
is obtained with variance 0{\) for each instance. Repetition can be used to reduce the variance 
below what is required. 

2. L-Ancilla Qubit Measurement Processes 

Often, we want to compute the expectation value of an operator O of the form 

M 



o = j2 a * u i v ^ w 



i=i 
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where Ui and Vi are unitary operators, a.j are real positive numbers (a,; > 0), and M is an integer 
power of 2. (In the case that M is less than a power of two, we can complete this definition 
by setting the a M+1 , ■ ■ ■ , a M > = 0, where M' is an integer power of two.) We can compute this 
expectation value by preparing M different circuits, each one with one ancilla qubit, and for each 
circuit measure (UjVi) (see Section |111 C 1[ ). Then, we multiply each result by the constant etj and 
sum the results. However, in most cases, the preparation of the initial state is very difficult. There 
is another way to measure this quantity by using only one circuit which reduces the difficulty. 
We first write the operator O as 

M 

0=MY,^U}V^ (27) 

i=i 

M M 

where N = ^2 a, and a? = a, /A/" (Yl a i = !)• Then we construct a quantum circuit with the 

t=i i=i 
following steps: 

1. Prepare the state |^ ) sucn that (^o\0^o) is the expectation value to be computed. 

2. Adjoin L ancillas to the initial state, where L = J + 1 and 2 J = M. The first of these 
ancillas, ai, is prepared in the state |+) = (|0) + |l))/\/2- This is done by applying the 
Hadamard gate to the initial state |0) (see Fig. f|a). The other ancillas, {a2, 23, • • ■ , a^} are 
kept in the state |0). 

3. Apply a unitary evolution E(a±, a 2 , • ■ ■ , «m) to the ancillas {a 2 , a 3 , • ■ ■ , a^} to obtain 

M 

\iP) = ai|00 - ■ - 0) + a 2 |00 - ■ -1) H ha M |ll • ■ ■ 1) = a { \i) , 

i=l 

where \i) is a tensorial product of the states (|0) or |1)) of each ancilla: \i) = \rj) a2 ® ■ ■ ■ ® 
\r]) 3L , where 77 can be or 1. The index i orders the orthonormal basis 

4. Apply the controlled unitary operations Ui which evolve the system by Ui if the state of the 
ancillas is |0) ai Then apply the controlled unitary operations Vi which evolve the system 
by Vi if the state of the ancillas is |l) ai |z). Once these evolution steps are finished, the state 
of the whole system is 

M M 



-7= \^)^Y.^\^)U l + \l) ai Y j a l \i)V l 

* t- i=l i=l 



l*0) • 
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5. Measure the expectation value of 2<7+ = o^ 1 + icr^ 1 = 2|0) ai (l|. It is easy to see that it 



M 



corresponds to the expectation value of the operator ^ of U]Vi. 



i=l 



6. Obtain the expectation value of O by multiplying (2cr+) by the constant M. 



The quantum circuit for this procedure is given in Fig. 10 



3. Measurement of Correlation Functions 

We now consider measuring correlation functions of the form Cab = (T^ATB), where T is a 
unitary operator and A and B are operators that are expressible as a sum of unitary operators: 

A = Y^ ouA and B = ^ (SjB 3 . (28) 

i 3 

The operator T is fixed by the type of correlation function that we want to evaluate. In the case of 
dynamical correlation functions, T is e~ tHt where H is the Hamiltonian of the system. For spatial 
correlation functions, T is the space translation operator e~ %p ' x (p and x are configuration space 
operators). The method for measuring these correlation functions is the same method described 
in Section pi C 1| or Section pi C 2[ We can use either the one- or the L-ancillas measurement 



process. 

To minimize the number of controlled operations and also the quantity of elementary gates 
involved, we choose LfJ = T^Ai and Vj = TBj. Now, we have to compute (UjVj). In Fig. [FT 



we show the circuit for measuring this quantity, where the circuit has only one ancilla in the state 
1+) = ( 1 ) + 1 1) ) / \/2- There, the controlled operations were reduced by noting that the operation 
of T controlled on the state |0) of the ancilla followed by the operation of T controlled on the state 
|1), results in a no-controlled T operation. This is a very useful algorithmic simplification. 



4. Measurement of the Spectrum of an Hermitian Operator 

Many times one is interested in determining the spectrum of an observable (Hermitian operator) 
Q, a particular case being the Hamiltonian H. Techniques for getting spectral information can be 
based on the quantum Fourier transform [ pPj , [H] ] and can be applied to physical problems [|2"2"|]. For 
our purposes, the methods of the previous Sections yield much simpler measurements without loss 
of spectral information. For a given H, the most common type of measurement is the computation 
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of its eigenvalues or at least its lowest eigenvalue (the ground state energy). To do this we start 
from an state \<p) that has a non-zero overlap with the eigenstates of H . (For example, if we want 
to compute the energy of the ground state, then \<p) has to have a non-zero overlap with the ground 
state.) For finite systems, \(j>) can be the solution of a mean-field theory (a Slater determinant in the 
case of fermions or Perelomov-Gilmore coherent states in the general case). Once we prepare this 
state (Section |lll A| and Appendix A), we compute (U{t)) = (<fi\U ~{t)4>) , where U is the evolution 
operator U (t) = e~ tHt . We then note that 

c 

|0> = 5>J*n>, (29) 

n=0 

with |\P n ) eigenstates of the Hamiltonian H. Consequently 

c 

(U(t)} = W 2 , (30) 

n=0 

where A n are the eigenvalues of H. The measurement of (U(t)) is easily done by the steps de- 
scribed in Section |111 C 1| (setting V = U(t) and U — 1 in Fig. Once we have this expectation 
value, we perform a classical fast Fourier transform (i.e., J (U(t))e tXt dt) and obtain the eigenvalues 
A n (see Appendix B): 

c 

FFT[(U(t))] = ^\ln\ 2 S(X - A re ) . (31) 

n=0 

Although we explained the method for the eigenvalues of H, the extension to any observable Q is 
straightforward, taking U (t) = e~ tQt and proceeding in the same way. 

Two comments are in order. The first refers to an algorithmic optimization and points to de- 
creasing the number of controlled operations (i.e., the number of elementary gates implemented). 
If we set V = e~ tQt , W — 1 (see Fig. ^) and perform the type of measurement described in 
Section pi C 1| the network has total evolution (ancilla plus system) e -l< ^e l< ^2, while if we set 
V — — e~ tQ z the total evolution is e* Qcr ^. Thus, this last algorithm reduces the number of 
gates by the number of gates it takes to represent the operator e~ lQ ^. The circuit is shown in 
Fig.0. 

The second comment refers to the complexity of the quantum algorithm as measured by system 
size. In general it is difficult to find a state whose overlap scales polynomially with system size. If 
one chooses a mean-field solution as the initial state, then the overlap decreases exponentially with 
the system size; this is a "signal problem" which also arises in probabilistic classical simulations of 
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quantum systems. The argument goes as follows: If \<f>) is a mean-field state for an iV d (=volume) 
system size whose (modulus of the) overlap with the true eigenstate is [7] < 1, and assuming that 
the typical correlation length of the problem £ is smaller than the linear dimension N, if we double 
N the new overlap is ~ e 2d ln ' 7 ' . 

We would like to mention that an alternative way of computing part of the spectrum of an 
Hermitian operator is using the adiabatic connection or Gell-Mann-Low theorem, an approach 
that has been described in [(]]] . 



5. Mixed and Exact Estimators 



We already explained how to compute different types of correlation functions. But in most 
cases, we do not know the state whose correlations we want to obtain. The most common case is 
wanting the correlations in the ground state \^ ) of some Hamiltonian H. Obtaining the ground 
state is a very difficult task; however, there are some useful methods to approximate these correla- 
tion functions. 

Suppose we are interested in the mean value of a unitary operator 0(t). If we can prepare the 
initial state \ ^f T ) in such a way that |^ ) = |^ / r) + e| < ^ ) ) (e is intended to be small), then after some 



algebraic manipulations [[23]], we have 

QPo|0(t)|tt T ) = 1 
(*o|*t> 2 



(#o|0(t)|tf > + (*T\0(t)\y T ) 



+ £>(e 2 ) (32) 



(*o|*o) (V t \Vt) 
where the term on the left-hand side of Eq. ^ is known as the "mixed estimator." Also, we can 
calculate the second term on the right-hand side of Eq. [32| with an efficient quantum algorithm, 
since we are able to prepare easily \^t)- Next, we show how to determine the mixed estimator 
using a quantum algorithm. 

If |\Po) is the ground state, then it is an eigenstate of the evolution operator U (£') = e~ lHt ' , and 
we can obtain the mixed estimator by measuring the mean value of U(t')0(t): Because \^t) = 
^2a n \^ n ) where a n = (^ n |^r) and |^ n ) are the eigenstates of H (U(t')\^ n ) = e~ iXnt '\^ n )) we 



can measure (Section IIIC3) 



(y T \U(t')0(t)\y T ) = Y,e lKt ' (^ T \^n)(^n\0(t)\^ T ) (33) 

n 

By performing a Fourier transform in the variable t' (F(u) = J e lu)t ' F{t')dt') in Eq. |33| and making 
the relation between the expectation value for time t and the expectation value for 0(t) = 11, we 
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obtain the value of the mixed estimator. Then, by using Eq. 1321, we obtain ™™ up to order 





By similar steps, we can obtain expectation values of the form m fig ,)" f° r a ^ n an ^ n '- The 




trick consists of measuring (Section [HI C 3| ) the mean value of the operator U(t')0(t)W(t") in 



the state \ <& T ) 



(y T \U(t')0(t)fr(t")\y T ) = J2e iXn *e iXn '*\*T\*n)(*n>\*T)(*n\0(t)\* n ,) (34) 



and then by performing a double Fourier transform in the variables t' and t" (F(A, A') = 
J e lXt ' e lXt F{t' ,t")dt'dt") we obtain the desired results. A particular case of this procedure is 
the direct computation of the exact estimator ^jjPfJ}^ ■ 



IV. APPLICATION TO FERMIONIC LATTICE SYSTEMS 

In this Section, we illustrate a procedure for simulating fermionic systems on a quantum com- 
puter, showing as a particular example how to obtain the energy spectrum of the Hubbard Hamil- 
tonian for a finite-sized system. We will obtain this spectrum through a simulation of a quantum 
computer on a classical computer, that is, by a quantum simulator. 

We start by noting that the spin-fermion connection described in Eq. [15] and Eq. [16| implies 
that the number of qubits involved in a two-dimensional lattice is L = 2(N X x N y ) if one uses 
the standard model to simulate spin-1/2 fermions. Also, the number of states for an L-qubits 
system is 2 L . From this mapping, the first N x x N y qubits represent the states which have spin-up 
fermions, and the other qubits ((N x X N y + 1) to 2(N X x N y )) spin-down fermions. In other 
words, if we have a system of 4 sites and have a state with one electron with spin up at the 
first site and one electron with spin down at the third site, then this state in second quantization is 
= a\.rCi\. \ | vac), where the fermionic operator aj. CT creates a fermion in the site j with spin a, 
and | vac) is the state with no particles (vacuum state). In the standard model, this state corresponds 



= ^n^+l^ = 1°) ® I 1 ) I 1 ) ® I 1 ) ® I 1 ) I 1 ) 1°) I 1 ) ITIUUU), (35) 



i=i 

where |vac) is the vacuum of the quantum spin 1/2, which we have chosen to be |||||||||). 
To represent the L-qubit system on a classical computer, we can build a one-to-one mapping 



to 
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between the 2 L possible states and the bit representation of an integer I defined by 

L 

I = £[n(i) x 2*- 1 ] (36) 
i=i 

where n(i) (occupancy) is if the spin of the i-qubit is |1) (J,), or 1 if the state is |0) (f). In this 
way, the state described in Eq. ^5] maps to I = 65. Because we are interested in obtaining some of 



the eigenvalues of the Hubbard model, we added an ancilla qubit (Fig. [12]). The "new" system has 
L = 2(N X x N y ) + 1 qubits, and we can perform the mapping in the same way described above. 

To simulate the evolution operator U(t) = e~ lHt on a classical computer using the above 
representation of quantum states, we programmed the "elementary" quantum gates of one and two 
qubits interactions. Each L-qubit state was represented by a linear combination of the integers I 
(Eq. |3(]). In this way, each unitary operation applied to one or two qubits modifies I by changing 
a bit. For example, if we flip the spin of the first qubit, the number / changes by 1 . 

We want to evaluate some eigenvalues of the spin- 1/2 Hubbard model in two spatial dimensions. 
The model is defined on a rectangle of N x x N y sites and is parametrized by spin preserving 
hoppings t x and t y between nearest neighbor sites, and an interaction hi on site between fermions 



of different ^-component of spin (Fig. |13J). The Hamiltonian is 

H = ~ t tx ( a \i,jy,<T a (i+l,3>,<r + a li+l,j);a a (i,J);o) + *w( a (ij);<T a (i,J+l);<7 + a \i,j+l);a a (i J);<r)] 

+ U n (i,3);l n (i,3);l ' ( 37 ) 

(i,3l 

where n^jy )(T = at .,. ff o^ -y a is the number operator, and the label (i, j);a identifies the site (X-Y 
position) and the ^-component of spin (er = ±1/2). We assume the fermionic operators satisfy 
strict periodic boundary conditions in both directions: aa,j);a = a {i+N x ,j);a = a (i,j+JV„);o-. 



To obtain the energy spectrum for this model, we use the method described in Section |IIIC 3 



(See Fig. |12| .) For this, we represent the system in the standard model, using the Jordan- Wigner 
transformation, mapping a two-dimensional spin- 1/2 system into a one-dimensional chain, with 
the use of Eq. [16] and Eq. [B| (Fig. §). 

As explained in Section pi C 4| , we find it convenient to start from the mean-field ground state 
solution of the model, represented by Hmf 

H M F = ~ [ txi^jy^i+lj);* + a (i+l,j);a a (i,j);<r) + *w( a (ij);«r a (ij+l);«T + a (i,j+l);<x a (i,j);a)] 

(id) 
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where the expressions in angular brackets are expectation values in the mean-field representation. 
Without loss of generality, we take U > and select the anti-ferromagnetic ground state mean- 
field solution. For this solution, we require N x and N y to be even numbers. If we were to simulate 
a one-dimensional lattice, we would however chose one of these numbers to be even and the other 
equal to 1. In the following we will only consider the half-filled case which corresponds to having 
one fermion per site; i.e., N e = N x x A^). 



First, we prepare the initial state. As discussed in Section pi A] , we do this by exploiting 
Thouless's theorem. We also use the first-order Trotter approximation (Section \ TB[ ), and then 
decompose each term of the evolution into one and two qubit interactions. Here, the matrix M 
now depends on the parameters of the Hamiltonian, as it does the ground state mean-field solution. 
After the decomposition, we then prepare the desired initial state by applying the unitary evolutions 
to a polarized state. (See Appendix A). 

Next, we execute the evolution U(t) = e~ lHt . For the sake of clarity we only present the 
first-order Trotter decomposition. To this end, we rewrote the Hubbard Hamiltonian as 

H = K + V = Ki + K^V , (38) 

where K a is the kinetic term (hopping elements with spin a) and V is the potential energy term. 
Because [./£„., V]_ ^ and = we approximated the short-time evolution operator 

U(At) by 

U(At) = e~ lHAt « e -iKM e -iVM. ( M ^ ) _ (39) 

N x xN y 

Because the term V = U Yl n (i,j);1 n (i,j);l = Vi 1S a sum °f operators local to each lattice 

site, each of these terms commute so 



e -iVAt 



l[e iVAt (40) 



The kinetic term is a sum over the bonds in the lattice (Fig. |13j): K a = K bond . a . Each bond 

bonds 



joins two nearest neighbor sites, either in the vertical or horizontal direction (Fig. [13]). Because 
of the periodic boundary conditions, the sites at the boundary of the lattice are also connected by 
bonds. We note that the terms in K that share a lattice site do not commute. For these terms we 
rewrite K a as 

K a = K°. a + K ia + K° y;a + K e ma , (41) 
where are the kinetic terms (for spin a) in the //-direction that involve the even (e) (and 



odd (o)) bonds in this direction (green and blue lines in Fig. [131). Then we perform the first-order 
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Trotter approximation 

e -iK a At ^ e -iK°. a At e -iKl^At e -iK°. a At e -iKl. a At _ ^ 

Because the odd and even bonds are not connected, each term in ( pdj ) is a sum of terms that 
commute with each other, that is: = K^ ,m , where [K^a' m , K e ^ ,m ]_ = 0, then: 

m 

e -iK<$At = JJ e -iK^ m At _ (43) 
m 

In summary we approximated the short-time evolution U (At) by 



e -iViAt 



(44) 



mi,m2,m3,m4;cr 

The total evolution operator is 

f>(t)=n^(At). (45) 

Each unitary factor in the evolution is easily decomposed into one and two qubit interactions 



(Section QB) 



The final step is the measurement process. To obtain some of the eigenvalues, we use the circuit 
described in Fig. [12|. Thus we are interested in the operator U(t/2)~°* instead of U(t/2) so we 
actually performed the first two steps after adding an ancilla qubit a (Fig. [12|), and then started with 
a "new" Hamiltonian H = —H <g> ^f, (and also a "new" evolution U(t) = e~ lHt ) and performed 
the same steps described above. 

The results for the simulation of the Hubbard model are shown in Fig. fRl. There, we also 



show the parameters At x and At 2 corresponding to the time-steps we used in the initial state 
preparation, where we used a first-order Trotter approximation, and in the time evolution, where 
we used a second-order Trotter approximation. 

In closing this Section we would like to emphasize that the simulation of the Hubbard model 
by a quantum computer which uses the standard model is just an example. Suppose one wants to 



simulate the Anderson model [ |24| ] instead using the same quantum computer, then similar steps 
to the ones described above should be followed. (There are two types of fermions but the iso- 
morphism still applies.) Similarly, if one wants to use a different quantum computer which has 
another natural "language" (i.e., a different operator algebra which therefore represents a different 
model of computation) one can still apply the ideas developed above simply by choosing the right 
isomorphism or "dictionary" |]3j]. 
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V. CONCLUDING REMARKS 



We addressed several broad issues associated with the simulation of physical phenomena by 
quantum networks. We first noted that in quantum mechanics the physical systems we want to 
simulate are expressed by operators satisfying certain algebras that may differ from the operators 
and the algebras associated with the physical system representing the quantum network used to 
do the simulation. We pointed out that rigorous mappings [|3|] between these two sets of operators 
exist and are sufficient to establish the equivalence of the different physical models to a universal 
model of quantum computation and the equivalence of different physical systems to that model. 

We also remarked that these mappings are insufficient for establishing that the quantum net- 
work can simulate any physical system efficiently even if the mappings between the systems only 
involves a polynomial number of steps. We argued that one must also demonstrate the main steps 
of initialization, evolution, and measurement all scale polynomially with complexity. More is 
needed than just having a large Hilbert space and inherent parallelism. Further, we noted that 
some types of measurements important to understanding physical phenomena lack effective quan- 
tum algorithms. 

In this paper we mainly explored various issues associated with efficient physical simulations 
by a quantum network, focusing on the construction of quantum network models for such com- 
putations. The main questions we addressed were how do we reduce the number of qubits and 
quantum gates needed for the simulation and how do we increase the amount of physical infor- 
mation measurable. We first summarized the quantum network representation of the standard 
model of quantum computation, discussing both one and multi-qubit circuits, and then recalled 
the connection between the spin and fermion representations. We next discussed the initializa- 
tion, evolution, and measurement processes. In each case we defined procedures simpler than the 
ones presented in our previous paper greatly improving the efficiency with which they can 
be done. We also gave algorithms that greatly expanded the types of measurements now possible. 
For example, besides certain correlation functions, the spectrum of operators, including the energy 
operator, is now possible. Our application of this technology to a system of lattice fermions and 
the construction of a simulator was also discussed and used the Hubbard model as an example. 
This application gave an explicit example of how the mapping between the operator of the phys- 
ical system of interest and those of the standard model of quantum computation work. We also 
gave details of how we implemented the initialization, evolution, and measurement steps of the 
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quantum network on a classical computer, thereby creating an quantum network simulator. 

Clearly, a number of challenges for the efficient simulation of physical systems on a quantum 
network remain. We are prioritizing our research on those issues associated with problems that 
are extremely difficult for quantum many-body scientists to solve on classical computers. There 
are no known efficient quantum algorithms for broad spectrum ground-state (zero temperature) 
and thermodynamics (finite temperature) measurements of correlations in quantum states. These 
measurements would help establish the phases of those states. Generating those states is itself a 
difficult task. 

Many problems in physics simulation, such as the challenging protein folding problem, are 
considered to be well modeled by classical physics. Can quantum networks be used to obtain 
significantly better (more efficient) algorithms for such essentially classical physics problems? 
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APPENDIX A: DIFFERENT STATE PREPARATION 



1. Coherent State Preparation: An example 



Here we illustrate by example the decomposition of an operator of the form e la Ma to generate 
an initial state. Typically M is generated by some mean-field solution to the physical problem of 
interest. Considerable detail is given. 

We consider 2 spinless fermions in a one-dimensional lattice of 4 sites (N e = 2 , n = 4). The 
operators a, and a] annihilate and create a fermion in the site j of the lattice. We want to prepare 



an initial state \<f)') = CqC^J vac) from the state 



OSr/2 



aja), | vac), where the operators Ck and c 



annihilate and create a fermion in the state of wave vector k, that is: 

3=1 



(Al) 



where k = 0, 7r/2, it, 3tv/2 are all possible wave vectors of the system, and Xj is the position in the 
lattice of the site (i.e., Xj = j — 1). 



From Eq. [AT], we see that the state \<f>') is a linear combination of states of the form a\a - |vac). 
The change of basis e lM (Eq. ^JJ) between the two sets of fermionic operators is: 



/ 4 \ 

Sr/2 



1 

2 



\ C 3tt/2 ' 



/I 
1 

1 

Vi 



1 \ 

—% 

-1 



(a\\ 



(A2) 



JMt 



If we calculate the eigenvalues and the eigenvectors of the matrix e lM , from Eq. [X2| we obtain 

f-1 0\ 
i 
10 

V o 001/ 

where Mr> is M in its diagonal form. Then, we have: 

/7T 



(A3) 



M, 



D 



-i log(e 



iMr 



0\ 

vr/2 



Vo o oo/ 



(A4) 
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To obtain the matrix M = A^M^A, we need to know the unitary matrix A, which is constructed 
with the eigenvectors of the matrix e tM . In this case we have: 



At 



/-1/2 l/y/2 1/V2\ 

1/2 -l/y/2 l/y/2 

1/2 -l/y/2 l/y/2 

\ 1/2 l/y/2 l/y/2 / 



(A5) 



hence, the Hermitian matrix M is: 





/ 1 


-1 


-1 




7T 


-1 


2 


1 





4 


-1 


1 


1 


1 




\-l 





1 


2 / 



In order to obtain we prepare the state |0) and then apply the evolution U = e 1 



(A6) 



If we 



want to simulate this fermionic system in a quantum computer (standard model), we have to use 
the spin-fermion connection (Section PQ ), and write the operator U as a combination of single 
qubit rotations and two qubit interactions. Also, the initial state \<j>) must be written in the standard 
model: 



a\a\ I vac) = <j\ 

|0)i® |0> 2 ®|1>3® 



^l(- ff M)l vac > = &l 



7^(j\ I vac) 



ITTU), 
|i)i ® |i) 2 



(A7) 
(A8) 

|1>« = IU---I) 



where the vacuum state in the standard model is |vac) 

i-i . _ 

((11 — o"i)o"i|vac) = aj|vac) = 0). With this mapping, the state \<f>') is a linear combination of 

1=1 

states of ^-component of spin 0. 



As noted in Section pi B| , sometimes the decomposition of the operator U in terms of one and 
two qubit operations is very difficult. To avoid this problem, we can use the Trotter decomposition 
(Eq. ^3]). In Fig. [15| we show the overlap (projection) between the state \<f>') and the state prepared 
using the first-order Trotter decomposition of U applied to the state \4>). 



2. Jastrow-type Wave Functions 



A Jastrow-type wave function is often a better approximation to the actual state of an inter- 
acting system, particularly when interactions are strong and short-ranged. Often one varies the 
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parameters in these functions to produce a state that satisfies a variational principle for some phys- 
ical quantity like the energy. Such states build in correlated many-body effects and are, in general, 
entangled states. The states described in the previous subsection (Appendix A 1) are unentangled. 
The classic form of a Jastrow-type wave function for fermions is[[T5]] 

|*o> = e s \4>') , (A9) 

where S = Yl a ij c l c j + Yl AjW c ! c ] c Jfc c J + • • • is an operator which creates particle and hole 

ij ijkl 

excitations, and \<p f ) is typically a Slater determinant. The iV-body correlations embodied in S 
take into account the short-range forces not included in | </>'). We will assume the and Pijki have 
been determined by some suitable means (for example, by a coupled-cluster calculation). If we 
decompose e s into a linear combination of unitary operators, we can then decompose |^o) mt0 
a linear combination of Slater determinants and thus prepare \ty ) as explained in [jlj]. Also, if 
the coefficients and are small, we can approximate e 5 by the first few terms in its Taylor 
expansion. Again, the state |^o) will be a linear combination of Slater determinants. 

Obviously, it is more natural for a quantum computer to generate a correlated state of the form 

|^o>=e*V>, (A10) 

where e lS is a unitary operator. In order to determine the iV-body correlation coefficients aty and 
Pijki, one could, in principle, use the technique of unitary transformations introduced by Villars 



26 



APPENDIX B: DISCRETE FOURIER TRANSFORMS 



In practice, to evaluate the discrete Fast Fourier Transform (DFFT) one uses discrete samples, 
therefore Eq. must be modified accordingly. In Fig. [13] we see that instead of having 5-functions 
(Dirac's functions), we have finite peaks in some range of energies, close to the eigenvalues of the 
Hamiltonian. Accordingly, one cannot determine the eigenvalues with the same accuracy as other 
numerical calculations. However, there are some methods that give the results more accurately 
from the DFFT. 

As a function of the frequency D, m , the DFFT (F(f2 m )) is given by: 



N-1 
3=0 



(Bl) 



where tj = jAt are the different times at which the function F is sampled (in the case of Section 



grgl F (tj) = (u(tj))), n r 



2nm 
NAt 



are the possible frequencies to evaluate the FFT of F(t) and 



N is the number of samples. (N must be an integer power of 2.) 

c 

Since we are interested in F(t) = Yl \ln\ 2 e~ lXnt (Eq. 



n=0 



N-1 



F(n m ) = At.J2hn\ 2 Ys el 



n=0 



and then 



3=0 



3 i(n m -A„)AiAT 
pi(Q m —\ n )At _ 



(B2) 



(B3) 



n=0 



If Q m is close to one of the eigenvalues A n and the A n are sufficiently far appart to be well resolved, 



we can neglect all terms in the sum other than n. If we take VL m and Vt 



m+l 



f2 m -\ 



2tt 
NAf 



both 



close to A n in such a way that \F(Vt m ) \ ^> \F(D, m+ i)\ 3> 0, then from Eq. |B3] we find that 



m+l, 



F(n r , 



gi(C m -A n )At _ j 
3 i(Q m+1 — X„)At 1 



After simple algebraic manipulations (and approximating ln(l + z) 
the correction to the energy A n : 

An = + AA n 



z for \z\ 



with 



AA„ w Re 



.F{Qm+l) 

1 F{n m ) 



D i2ir/N _ Y 

At 



(B4) 
0) we obtain 

(B5) 
(B6) 
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FIG. 1 : Relationship between different models of computation (with their associated operator algebras) and 
different physical systems. Question marks refer to the present lack of a quantum computer device using the 
corresponding elementary physical components indicated in the box. Diamond-shaped arrows represent the 
natural connection between physical system and operator language, while arrows on the circle indicate the 
existence of isomorphisms of *-algebras, therefore, the corresponding simulation of one physical system by 
another. A wave function view of this relationship is given in [0J. 
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FIG. 2: Bloch-Sphere representation of a one qubit state parametrized as |a) = cos ||0) + e lv sin The 
curved arrows indicate the sign of rotation of = R^—t) about the particular axis fj,. Our (arrow) 

color convention is: |0) — ► blue; |1) — > red; other linear combinations — ► magenta. 
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FIG. 3: (a) Some one qubit elementary gates and (b) a two qubit elementary gate. 
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(a) 




FIG. 4: (a) Hadamard gate decomposition and (b) Bloch-Sphere representation of a Hadamard gate applied 
to the state |+). 
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FIG. 5 : C-NOT gate decomposition and its matrix representation. The control qubit is 1 . Note that the last 
circuit realizes the C-NOT matrix operation up to a global phase e~ 1 ^. 
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FIG. 6: Bloch-Sphere representation of the state obtained by the C-NOT gate applied to the "classical" state 
1 10). The sequence of elementary operations is the same as fig. |5| (time flows from left to right with the 
lower row continuing the upper one). For each Bloch-Sphere the two arrows indicate the states of the two 
qubits, with the left representing qubit one. 
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FIG. 7: (a) C-U operation with the state of the control qubit a being in | l) a and (b) C-U' operation controlled 
with the state |0) a . (See text for notation.) 
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FIG. 8: Mapping used to connect the labels of a two-dimensional N x x N y lattice to the labels of a chain 
(i.e., a one-dimensional array of integer numbers). 
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FIG. 9: Measurement of physical quantities using one extra (ancilla) qubit |a). In this case (2a* 
(* \rtV\* ). 
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FIG. 10: Measurement of physical quantities using L-ancillas qubits {ai, ■ • • , ax}. In this case (cl 1 ) = 

M 

2^(*o|[ Eoi^Ml l*o) (see text). 
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FIG. 11: Circuit for the measurement of spatial and time correlation functions. In this case {2a a + ) = 
(T* AiTBj). Notice the simplification achieved by reducing two C-T operations into only one uncontrolled 
T operation. 
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FIG. 12: Circuit for the measurement of the spectrum of an Hermitian operator Q. In this case (2o+) = 
((/>|e"^V} (see text). 



41 




FIG. 13: Two-dimensional lattice in the Hubbard model. Here, the green and blue arrows identify the even 
and odd bonds. 
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FIG. 14: Energy spectrum of the Hubbard model obtained from the simulator. The lattice has 4x2 sites 
(which requires 16 qubits), with t x = 1, t y = 1 and U = 4 and the time steps used in the Trotter approxi- 
mation (to prepare the initial state and apply the evolution) are At± = At2 = 0.05. 
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FIG. 15: Overlap between the exact initial state and the state prepared with the Trotter decomposition, for a 
system with two fermions in a 4-sites lattice. 
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